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Abstract 


Rheological parameters are the key factor for hydraulic calculations and hole cleaning efficiency in the drilling 
operations. Rheological parameters include mud density, plastic viscosity, apparent viscosity, and yield point. Rheological 
parameters should be measured or estimated during the real operation to determine the actual equivalent mud 
circulation density, surge and swab pressures, hydraulic calculations, and hole cleaning efficiency. The mud rheological 
properties can be obtained from laboratory measurements. These measurement scan be divided into two types; the first 
one is simple testing such as mud density and Marsh funnel viscosity. The second set of experiments is the rheological 
properties measurements. Onsite, Marsh funnel viscosity and fluid density are measured frequently (every 10 to 20 
minutes) to monitor the change of the drilling fluid properties while as the rheological properties are measured twice a 
day. 

In this paper a new model was developed to estimate the rheological parameters of invert emulsion-based drilling fluid 
using Marsh funnel viscosity and mud density. Multiple nonlinear regression technique was used to develop the new 
model. 6,000 data points from field measurements of the Marsh funnel viscosity, mud density, plastic viscosity, and yield 
point, were used in this study to develop the new model. The new model can be used to predict the following properties 
of the drilling fluid; viscometer readings at 300 and 600, plastic viscosity, yield point, apparent viscosity, power law 
index, and consistency index. All the previous parameters will be predicted using only Marsh funnel time and drilling 
fluid density. The developed models consists of six correlations that can be used to estimate the rheological properties 
with an average absolute error less than 4% and a correlation coefficient greater than 0.95 based on the Marsh funnel 
time and drilling fluid density. These correlations can be used to prepare the input parameters for hole cleaning and 
hydraulic calculations in real time. Applying this inexpensive technique will help drilling engineers predict any problem 


before it happens such as pipe sticking, sage issues, and loss of circulation. 
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Introduction 


P Drilling fluid is the key element in the drilling 
operations. It has many functions such as controlling the 
formation pressure, carrying the drilled-cutting form 
down hole to the surface, transferring information from 
down hole to surface, cooling andlubricating the drill 
string and the drill bit, and form ideal filter cake for 
wellbore stability purposes. In order to perform these 
functions, the drilling fluid should contain inert solids for 
density control, reactive solids for viscosity control, and 
polymers for filtration control [1]. This complex non- 
Newtonian fluid mixture will be used to drill different 
wells, such as vertical, deviated, horizontal, and extended 
reach wells. 


Invert emulsion is a type of oil based drilling fluid that 
has 50% water ratio. This type of oil is low toxicity and oil 
is the continuous phase while water with CaCl2 is the 
dispersed phase. Calcium chloride is used to increase the 
salinity of water and prevent any migration of water 
molecules from solution to formation. Invert emulsion is 
commonly used in the drilling operations in oil and gas 
industry [2]. 


Different instruments are used to determine the fluid 
properties, such as mud balance, Fann 35 viscometer, API 
filter press, Marsh Funnel, and high-pressure high 
temperature filter press. These measurements can be 
performed in the laboratory and it takes long time, on the 
well site, some of these instruments are not available and 
will take long time which is not available during the real 
drilling operation. The drilling fluid should be designed to 
have viscosity to sweep the cutting to the surface. This 
term of viscosity is called plastic viscosity (PV) according 
to Bingham model and consistency index (k) according to 
Power low model. In addition, the drilling fluid should 
have a yield stress which is called yield point (YP) 
according to Bingham model and flow behavior index (n) 
according to Power Low model.Drilling fluid rheology 
plays a key factor in many aspects of the drilling 
operations. It can affect the following parameters; 
hydraulic calculations, hole cleaning, and sag tendency of 
weighting material (separation and settling of the 
weighting material), Tehrani [3]. The ratio of YP/PV is 
very important for hole cleaning and preventing common 
problems in the drilling operations such as pipe sticking, 
loss of circulation, and sagging effect. The drilling fluid 
should be designed to prevent any sag issue and fluid loss 
problem [4,5]. The formation damage can be minimized 
by using microemuslion; oil based drilling fluids, and 
inverted emulsion drilling fluid [6]. The invasion of the 
inverted emulsion drilling fluid system is very low 
compared to the water based drilling fluids because of the 
high viscosity of the oil and this will minimize the damage 
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caused by the drilling fluid invasion. In turn, the need for 
well stimulation operation after drilling will be minimized 
if not eliminated. 


The knowledge of rheological properties of the drilling 
fluid, which are generated from suspended solids and 
rheology control additives, is highly important in the 
drilling operations [7]. Each parameter plays an 
important role in the success of drilling operations. 
Zamora and Power stated that the yield stress is the most 
important factor for hole cleaning problems, sag issue, 
and surge and swab pressures [8]. They introduced a 
relation between the yield stress (ty) and the yield point 
(YP) for different drilling fluid types. They concluded that 
the acceptable range of ty/YP ratio is 0.5 - 0.68 for 
synthetic-based mud, 0.48 - 0.59 for oil-based mud, and 
0.2 - 0.4 for water-based mud. The t,y/YP is an important 
parameter to describe the mud rheology.The normal 
operation on the well site is to perform a complete test for 
the drilling fluid at morning and evening. Only mud 
density and Marsh funnel viscosity are measured 
frequently to give an indication of the change in drilling 
fluid properties. Marsh funnel is used to obtain the time in 
seconds required to fill a 950 cm? volume cup with 
drilling fluid [9]. Marsh funnel can be used to estimate 
many parameters such as fluid yield stress. This 
technique is inexpensive and it gives accurate results as 
Fann 35, [10]. The determination of rheological 
parameters is very important for pressure loss 
calculations. Other factors should be taken into 
consideration for hydraulic calculations such as drill 
cuttings. Zhang studied the effect of solid content in the 
calculation of the pressure loss in the annulus. Zhang et al. 
concluded that at low flow rates, the drilling fluid solid 
contentaffected the pressure profile while at high flow 
rateitseffect diminished [11]. Due to this effect, the 
bottom hole pressure may increase when the pump 
stopped. Solid effect should not be neglected in hydraulic 
calculations because it maycause pressure loss greater 
than the friction pressure loss. In addition to that, the hole 
shape has high effecton hydraulic performance and hole 
cleaning. Noncircular holes may have lower pressure 
gradient and better hole cleaning than circular holes [12]. 
In the literature, only few correlations were found to 
determine the apparent viscosity from the rheometer or 
Marsh funnel viscosity. In this paper we used 6000 data 
points collected from actual field measurements of invert 
emulsion drilling fluid to develop the correlations to 
determine the drilling fluid rheological properties. Also 
we will compare the developed correlations with the ones 
existed in the literature to show the accuracy of the 
developed correlations. 
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Apparent Viscosity from Marsh Funnel 


Based on 26 drilling fluids, a graphical method was 
developed to estimate the apparent viscosity of these 
drilling fluids as a function of Marsh funnel and the 
Stormer viscometer [13]. Equation 1 shows the 
relationship of the apparent viscosity as a function of 
Stormer speed and some constants. 


K 
AV == (1) 
AV apparent viscosity, cP 
R Speed of Stormer cylinder, rpm 
x exponent related to initial shear stress 


constant (x and K can be determined from 
the plot of Stormer viscosity measurement vs the 
viscometer speed) 


Pitt developed a correlation to convert Marsh funnel 
time to effective viscosity of drilling fluids as a function of 
drainage time and mud weight. He built a numerical 
simulation model of the flow behavior of the drilling fluid 
in Marsh funnel was built. The following correlation was 
developed to estimate the effective viscosity of the drilling 
fluid, Equation 2 [14]. 


AV = D(T — 25) (2) 


Almahdawi introduced a new model to estimate the 
apparent viscosity of drilling fluids using Marsh funnel 
time and fluid density [15]. Several samples of drilling 
fluids were prepared with different quantities of barite, 
bentonite, and other additives in order to give various 
rheological behaviors. Fluid density and Marsh funnel 
time of the samples were measured. The results were 
used to fit Equations 3 and 4. They concluded that 
Equation 3 related the apparent viscosity to Marsh funnel 
time gives has accurate results compared to Equation 4. 
They found out that the constant of 28 is more 
appropriate than 25, obtained by Pitt [14]. 

AV = —0.0118 « T? + 1.6175 * T — 32.168 (3) 
AV = DCT — 28) (4) 


Where 
AV apparent viscosity, cP 
D Fluid density, g/cm? 
T Marsh funnel time, sec. 


The objective of this paper is to develop new empirical 
correlations of plastic viscosity, yield point, apparent 
viscosity, flow behavior index, and flow consistency index. 
These correlations will be developed for invert emulsion- 
based drilling fluid based on the analysis of 6000 actual 


Elkatatny SM, et al. Determination of the Rheological Parameters of Inve 


field data measurements. Table 1 lists the main available 
parameters for this study. 


Density T PV YP 

: (lb/ft?) |(second/quart)| (cP) | (lb/100 ft?) 
1 57 44 11 16 
2 57 45 11 19 
7 70 65 22 21 
8 71 70 26 23 
9 72 71 22 21 
6000} 110 85 40 26 


Table 1: Field data for invert emulsion based mud. 


The importance of this research is the real time 
determination (every 10 to 15 minutes)of the hydraulic 
calculations, surge and swab pressure, hole cleaning 
efficiency, and equivalent circulation density based on the 
rheological parameters that will be determined from the 
mud density and Marsh funnel viscosity. Real time 
monitoring of different drilling parameters will help early 
detect and mitigate the possible drilling problems. 


Nonlinear Curve Fitting 


A complete analysis for the invert emulsion drilling 
fluid was performed to predict the rheological properties 
of the drilling fluid and we compared the estimated values 
to the given value. Apparent viscosity, the viscometer 
reading at 600 rpm (Reoo), and the viscometer reading at 
300 rpm (R300) can be determined from the available data 
using Equations 5 to 7. 


R300 = YP + PV (5) 
Reoo = PV + Roo (6) 


R 
AV = 3°) 


The flow behavior index and flow consistency index are 
calculated using Equations 8 and 9. 


R 
n = 3.32 * log (=) (8) 
300 
(9) 


_ _‘*eo0 
1022” 


Using non-linear curve fitting technique, Equation10 
was developed to predict the viscometer reading at 300 
rpm using the drilling fluid density and Marsh funnel 
viscosity.The average absolute error of the estimated 
values of R300 was calculated based on the actual one and 
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it was 2.82 and the correlation coefficient was 0.96 (R2 = 7] 
0.92), as shown in Figure 1. 
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20 Actual Values of Fluid Behavior Index 
20 40 60 80 Figure 3: Estimated value of flow behavior index vs. actual 
Actual Values of R309 value using nonlinear regression technique. 
Figure 1: Estimated value of R300 vs. actual value of R300 
using nonlinear regression technique. R399 = 19.0001 * D * log(T — 35.748) + 12.896 (10) 
Roo0 = 37.2748 * D * log(T — 30.80) + 1.1638 (11) 
The same technique was used to estimate the n = 0.1193 * D * log(T — 11.249) + 0.3459 (12) 
viscometer reading at 600 rpm based on the fluid density = 
and Marsh funnel time. Equation 11 was developed to D drilling fluid density, g/cm? 
estimate Reoo with an average absolute error of 3.28 anda T Marsh Funnel time, second/quart 
correlation coefficient of 0.96 (R2 = 0.93), Figure 2. In 
order to predict the parameters for power low model, Validation of Developed Correlations 
Equation 12 was developed based on fluid density and 
Marsh funnel time. To validate the developed correlations, plastic viscosity 


and flow consistency index were calculated based on the 
estimated values of R300, Reoo, and flow behavior index. 
The obtained results were plotted vs. the actual measured 
heatian Machin Oar a8 data to investigate the accuracy of the developed 
; TR? =0.93) correlations. Plastic viscosity was calculated based on 
Equation 6. Figure 4a shows the comparison of the 
calculated values of plastic viscosity based on estimated 
values of R300 and Reoowith the actual values of plastic 
viscosity. 


Estimated Values of Reoo 


[ Average Absolute Error = 4.0 


20 40 60 80 100 120 140 R? = 0.94 


Actual Values of Ryoo 


Figure 2: Estimated value of Reoo vs. actual value of Reoo 
using nonlinear regression technique. 


Figure 3 shows the comparison of estimated values of n 
with the actual ones with an average absolute error of 1.4 
and a correlation coefficient of 0.96 (R2 = 0.92). 20 30 40 

Actual Values of Plastic Viscosity, cP 
Plastic viscosity from estimated values of 
viscometer reading. 


Calculated Values of Plastic Viscosity, cP 
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Average Absolute Error = 3.68 


R? =0.90 


Estimated Values of Consistency Index 


0.3 0.5 0.7 0.9 + 13 15 
Actual Values of Consistency Index 


b. Consistency index from estimated values of 
Reo and flow behavior index. 


Figure 4: Validation of the developed correlation by 
calculating the plastic viscosity and consistency index. 


It was found that the average absolute error was 4.0% 
and the correlation coefficient was 0.97 (R2 = 0.94) 
comparing the actual and calculated values of the plastic 
viscosity. Flow consistency index was calculated based on 
the estimated value of flow behavior index and 
viscometer reading at 600 rpm, Equation 9. Figure 4b 
shows that the average absolute error of 3.68%and the 


correlation coefficient of 0.95 (R2= 0.90) comparing the 
calculated value of flow consistency index with the actual 
one. 


To estimate the yield point, plastic viscosity, and 
apparent viscosity directly from the Marsh funnel time 
and fluid density, non-linear regression technique was 
used and to obtain Equations 13-15 for yield point, plastic 
viscosity, and apparent viscosity, respectively. The same 
results were obtained using Equations 5-7 based on the 
estimated values of the rheometer readings. The 
importance of developed equations (13-15) is the 
directed estimation of these parameters from the Marsh 
funnel viscosity and mud density which will make the 
programming of these models easy and simple for the real 
time prediction. 


YP = 1.5492 « D = log(T — 30.9) + 18.84 (13) 
PV = 17.735 * D * log(T — 30.9) — 7.979 (14) 
AV = 18.833 « D + log(T — 30.9) + 0.9186 (15) 


D drilling fluid density, g/cm3 

T Marsh Funnel time, second/quart 
YP yield point, lb/100 ft2 

PV plastic viscosity, cP 
AV apparent viscosity, cP 


Fluid Density, lb/ft? 


Fluid Density, lb/ft’ 


Yield Point, Ib/100 ft? 
as ” ” 2 = os ” oo o2 


8 3 


8 


Fluld Density, lb/ft’ 
Fluid Density, tb/ft* 
3 


$ 


Flow Behavior index 


Plastic Viscosity, cP 
wo 40 


~” 20 


Fluid Density, lb/ft’ 
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Figure 5: Comparison of the rheological parameters over a wide range of fluid density. 
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Figure 5 shows the comparison of the predicted values 
of the rheological parameters with the actual ones for a 
wide range of density (60 to 110 lb/ft). Once the 
rheological parameters are available, Bingham model or 
Power low model can be usedfor drillinghydraulics 
calculations. Knowing the actual drilling parameters, such 
as rate of penetration, the pipe rotational speed, flow rate, 
and pump pressure, the hole cleaning efficiency can be 
determined based on the estimated values of the 
rheological parameters. 


Validating the Apparent Viscosity Correlations 
with Previous Correlations 


As explained in the introduction part, there are only 
two available correlations in the literature that discussed 
the prediction of the apparent viscosity from the mud 
density and Marsh funnel viscosity. In this section we will 
compare the developed correlation of the apparent 
viscosity with Pitt’s and Almahdawi et al. correlations 
[14,15]. 


& 


$ 


Estimated Values of Apparent Viscosity, cP 


Actual Values of Apparent Viscosity, cP 


a. Apparent viscosity by the new developed 


correlation 


Average Absolute Error = 13.6 
Ro 


) 
i 


60 80 100 
Actual Valves of Apparent Viscosity, cP 


b. Apparent viscosity 


Figure 6: Comparison of estimated values of apparent 
viscosity by the new correlation with model. 
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Figure 6 shows that the average absolute error was 
3.28% and the correlation coefficient was 0.96 (R? = 0.3) 
using Equation15 to estimate the apparent viscosity 
based on mud density and Marsh funnel viscosity. Pitt 
correlation overestimated the apparent viscosity as the 
average absolute error was 13.6% and the correlation 
coefficient was 0.84 (R2 = 0.71). Almahdawi et al. 
correlation yielded an average absolute error of 43%. 
Based on these results, it can be concluded that the 
developed correlation is accurate enough to predict the 
apparent viscosity with the highest correlation coefficient 
and the lowest average absolute error compared to 
previously developed correlation. This is the first time to 
predict other rheological parameters such as plastic 
viscosity, yield point, flow behavior index, and flow 
consistency index based on mud balance and Marsh 
funnel viscosity. The developed correlations have a high 
correlation coefficient (> 0.95) and low average absolute 
error (<4%). 


The Importance of Predicting Drilling Fluid 
Properties in Real Time 


The real time prediction of power law index, 
consistency index, plastic viscosity, yield point, and 
apparent viscosity from the Marsh funnel and drilling 
fluid density, will help monitor the drilling process. For 
example the hole cleaning efficiency can be determined 
during drilling in real time based on the predicted 
parameters using the developed model. The prediction of 
hole cleaning efficiency will prevent several drilling 
problems such as high drag, torque, and pipe sticking. 
Once the power law index is determined and consistency 
index in real time using the developed model, the cutting 
carrying efficiency can be determined as follows [16]: 


a 
n = 
Mine = ‘me 


[sezexto" re 1) (a (2) ie (18) 
Pmud Dp—Dp 


Where; yn = hole cleaning efficiency; V. = drilling fluid 
critical velocity; Vme = drilling fluid effective velocity; n = 
power law index; Pmua = drilling fluid density, Dp = hole 
size; Dp = drill pipe size. 


2n+1 


7 47) 


Using the developed model will add more dimensions 
to predict several factors that control the drilling 
operations such as hole cleaning efficiency, cutting 
concentration in the hole, and prediction of drilling 
problems. There is no current model that can predict the 
drilling parameters in real time (every 10 to 20 min) 
because all models depend on _ the laboratory 
measurements which are performed maximum twice a 
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day. Using the developed model we can predict real time 
and this can be logged with the well depth, any anomalies 
can be detected and mitigated immediately. This will save 
the cost of drilling operations and also will minimize the 
drilling time. 


Conclusion 


Non-linear regression technique was used to develop 
six correlations to determine the rheological parameters 
based on 6000 actual field measurements. The obtained 
correlations can be used to estimate the rheological 
parameters based only on Marsh funnel time and drilling 
fluid density with an average absolute error less than 4% 
and a correlation coefficient greater than 0.95. Drilling 
hydraulics and hole cleaning efficiency can be obtained 
and predicted in real time based on the developed 
correlations. This technique is inexpensive and will save 
time and will reduce the total cost of the drilling 
operations. 
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